Astronomy <fc Astrophysics manuscript no. 1350 


February 2, 2008 


(DOI: will be inserted by hand later) 





NGC 4254: a spiral galaxy entering the Virgo cluster 

B. Vollmer 1 ' 2 , W. Huchtmeier 2 & W. van Driel 3 

1 CDS, Observatoire astronomique de Strasbourg, UMR 7550, 11, rue de l'universite, 67000 Strasbourg, France 

2 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 

3 Observatoire de Paris, Section de Meudon, GEPI, CNRS UMR 8111 and Universite de Paris 7, 5 place Jules 
Janssen, 92195 Meudon, Cedex, France. 

Received / Accepted 

Abstract. Deep Effelsberg Hi spectra of the one-armed, bright Virgo cluster spiral galaxy NGC 4254 are presented. 
Five different positions were observed in the 21 cm Hi line with the Effelsberg 100-m telescope: one in the center 
and 4 located one HPBW to the NE, NW, SW, and SE, respectively, from the galaxy center. The spectra 
are compared to existing deep VLA observations, and the single dish and interferometric Hi data are used to 
constrain a dynamical model which includes the effects of ram pressure. The peculiar, one-armed spiral pattern 
of NGC 4254 and its distorted and kinematically perturbed atomic gas distribution can be explained by a close 
and rapid encounter ~ 280 Myr ago with another massive Virgo galaxy, followed by ram pressure stripping that 
is still ongoing. The stripping occurs almost face-on, since the angle between the disk and the orbital plane is 70° . 
The galaxy with which NGC 4254 had its encounter is tentatively identified as the lenticular NGC 4262. 
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dynamics 



1. Introduction 

NGC 4254 is a bright spiral galaxy which is located at the 
periphery of the Virgo cluster, at a projected distance of 
3.7° (~ 1 Mpc 1 ) from the cluster center (M87). Its basic 
parameters are listed in Tab. ^ The spiral structure of 
NGC 4254 shows an important m=l mode, giving it a one- 
armed appearance (Fig. [EJ) . The m=2 mode, which is most 
prominent in most grand-design spirals, is almost absent. 
Iye et al. (1982) suggested that the intrinsic properties 
of NGC 4254 are responsible for the predominant m=l 
mode. 

Phookun et al. (1993) obtained deep Hi 21 cm line ob- 
servations of this galaxy with the VLA. These deep obser- 
vations revealed, in addition to the galactic Hi disk com- 
ponent, Hi clouds superposed on and beyond the Hi disk 
with velocities that do not follow the disk rotation pattern. 
They also found very low surface density atomic gas up 
to ~ 6' (~ 30 kpc) to the north-west of the galaxy center. 
Based on their observations, Phookun et al. (1993) sug- 
gested that the low surface density gas represents a rem- 
nant of an entity that was tidally disrupted by NGC 4254 
and is now merging with it. In this scenario the m=l spiral 
structure was stimulated by the infalling gas. 

Send offprint requests to: B. Vollmer, e-mail: 
bvollmer@astro.u-strasbg.fr 
1 We use D = 17 Mpc for the distance to the Virgo cluster. 



Table 1. Physical Parameters of NGC 4254 



Other names M 99 

VCC 307 
CGCG 098 - 144 

a (2000) a 12 h 18 m 49.6 s 

S (2000) a 14°24'59" 

Morphological type a Sc 

Distance to the cluster center (°) 3.7 

Optical diameter D 25 a (') 5.4 

~ 10.44 

Systemic heliocentric velocity a (kms -1 ) 2407±3 

Distance D (Mpc) 17 

Vrot max (kms" 1 ) 150 b 

Majoraxisp.a. 68° b 

Inclination angle 42° b 

HI deficiency" 0.17±0.2 



a RC3, de Vaucouleurs et al. (1991) 
b Phookun et al. (1993) 
c Cayatte et al. (1994) 



NGC 4254 resembles in many respects NGC 4654 
(Phookun & Mundy 1995), another bright Virgo cluster 
spiral galaxy that is located at about the same projected 
distance from the cluster center, but at the opposite side 
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Table 2. Integration times and rms. resulting noise (Table |2J), which is partly determined by 

small amplitude interferences, is close to the theoretical 

noise of 2 mJy per hour of integration. 

position Center Northwest Southwest Southeast Northeast 

At (min) 60 60 60 40 45 3. Results 



rms (mJy) 2.8 1.8 1.4 3.5 1.7 



The 5 Effelsberg spectra are shown in Fig. [2] and Fig. 1111 
In order to compare our single-dish Hi spectra to inter- 
ferometric data, where the galaxy is spatially resolved, we 
use the VLA 21 cm C and D array data of Phookun et al. 
(1993). These data have a spatial resolution of 32" x 28" 
and a channel separation of 10 kms -1 . The bulk of the 
Hi emission detected with the VLA is contained within a 
~ 6' x 5' area, while a faint blob of Hi was found ~ 3' 
south-east of the galaxy center and a very faint Hi tail 
was detected towards the north-west, extending out to 10' 
from the center. 

In order to enhance the diffuse and extended Hi emis- 
sion, we smoothed the data cube to a resolution of 2', 3', 4', 
and 5'. Three dimensional masks were produced by setting 
all pixels to one where the intensity is above 5er and zero 
otherwise. This mask was applied to the initial data cube. 
For the 5 positions synthesized spectra from the VLA data 
were produced using the position of the pointing and the 
beamsize of the Effelsberg 100-m telescope at 21 cm (9.3'). 
We thus use the spatial information of the VLA data to 
construct the synthesized spectra. It turned out that a 
smoothing to a 4' resolution leads to a maximum flux den- 
sity in the spectra. Therefore, all synthesized spectra are 
made with a 4' smoothing. Moment maps were made of 
the masked cube without clipping (Fig.HJ. In these maps 
the north-western, very low surface density plume is well 
detected. 

The resulting spectrum for the central position is rep- 
resented by the dashed line in the lower panel of Fig. El 
We find a total single-dish flux of 67.7 Jykms -1 , which 
corresponds to an atomic gas mass of 4.6 10 9 M Q . The 
central velocity is 2406±10 kms -1 , the linewidths at 50% 
and 20% of the peak flux are W50 = 230 kms -1 and 
W20 = 270 kms -1 . The total and peak flux densities 
of Phookun et al. (1993) are within ~10% of our val- 
ues, which is within the calibration accuracy. Both spectra 
agree reasonably well. 

Fig-Elshows the 100-m spectra of the 4 off-center posi- 
tions (solid lines) together with the synthesized VLA spec- 
tra. The NE and SE Effelsberg and VLA spectra are in 
very good agreement. The Effelsberg NW spectrum shows 
excess emission between —100 and kms -1 and less emis- 
sion between and 100 kms -1 with respect to the VLA 
spectrum. The Effelsberg SW spectrum shows excess emis- 
sion around —200 kms -1 and less emission between and 
100 kms -1 with respect to the VLA spectrum. These dif- 
ferences can be explained by three effects: (i) pointing 
errors of the Effelsberg telescope, (ii) disk emission de- 
tected by the sidelobes of the Effelsberg telescope and 
(iii) extended (> 15') and/or very low surface density 
(< 2 x 10 18 cm -2 ) gas which has not been detected by 
the VLA. The excess emission of the NW and SW VLA 



of the cluster, since both galaxies (i) show a prominent 
m=l spiral structure, (ii) are not Hi deficient, and (iii) 
have an extended, very low surface density gas tail. 

Vollmer (2003) showed that the m=l spiral structure 
of NGC 4654 can be explained by a tidal interaction of 
NGC 4654 with its companion NGC 4639 that took place 
~ 500 Myr ago, while its extended gas tail can be ex- 
plained by ongoing ram pressure stripping (Gunn & Gott 
1972). Since the outer part of the gas disk was perturbed 
by the gravitational interaction, only a small amount of 
ram pressure is needed to create the gas tail. Thus, only 
the combination of a close encounter with NGC 4639 fol- 
lowed by continuing moderate ram pressure can reproduce 
all observed properties of NGC 4654. In this article we in- 
vestigate if the same evolutionary scenario, involving a 
close encounter and ongoing ram pressure stripping ap- 
plies to NGC 4254. 

The Effelsberg observations are described in Sect. [21 
followed by the results in Sect. [31 The model is introduced 
in Sect. 01 We present two different simulations in Sect. El 
(i) a close and rapid encounter between NGC 4254 and 
another massive Virgo cluster galaxy (Sect. I5.1fl and (ii) 
a close and rapid encounter together with ram pressure 
stripping (Sect. l5~2*|l . The simulations are compared to Hi 
observations in Sect. 15.31 and discussed in Sect. El followed 
by our conclusions in Sect.0 

2. Observations 

Between January and March 2002 we performed 21-cm 
line observations with the 100m Effelsberg telescope at 5 
different positions in and around NGC 4254 at a frequency 
corresponding to the systemic velocity of NGC 4254 and 
a bandwidth of 12.5 MHz. The two-channel receiver had 
a system noise of ~30 K. The 1024 channel autocor- 
relator was split into four banks of 256 channels each, 
yielding a channel separation of ~10 kms -1 . We further 
binned the channels to obtain a final channel separation 
of ~20 km s _1 . The central position of the galaxy and four 
positions at a distance of one HPBW (9.3') to the NW, 
SW, SE, and NE, respectively, from the galaxy center were 
observed in on-off mode using pairs of 5 min on source and 
5 min off source integrations. The total integration times 
are shown in Tabled Care was taken to avoid other Virgo 
galaxies with radial velocities within our bandwidth in all 
5 on and off source positions. We used 3C286 for point- 
ing and flux calibration. The on and off observation times 
ranged between 80 min and 120 min per position. The 
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Fig. 2. Solid lines: Effelsberg 100-m spectra of the four off-center positions. Their locations with respect to the galaxy 
center are marked on top of each panel. Dashed line: synthesized VLA spectra which only show Hi disk emission. 
Dotted line: 3cr noise levels of the Effelsberg spectra. Velocities are given relative to the systemic velocity of NGC 4254, 
2407 kins" 1 . 



spectra in the velocity range between and 100 km s^ 1 are 
possibly due to pointing errors of the Effelsberg telescope. 
A pointing error of about 10% of the telescope's HPBW 
to the west can account for the VLA excess emission with 
respect to the Effelsberg spectrum. However, such a point- 
ing error is unlikely, though in principle not impossible, as 
it is 3 to 5 times the typical pointing correction we had to 
apply during the observing run. It is not excluded that the 
Effelsberg excess emission with respect to the VLA spec- 
trum is due to a sidelobe detection of the disk emission of 
the galaxy. The Effelsberg antenna diagram shows 4 first 
sidelobe maxima at a distance of ~ 13' from the beam cen- 
ter and at a -20 db level (1 % of the peak). One percent 
of the emission of the central profile (Fig. 1111) is ~ 3 mJy, 
which corresponds approximately to the 3ct noise levels 
indicated in the NW and SW spectra. Indeed, an inspec- 
tion of the parallactic angle of the Effelsberg observations 
showed that its variation was small for the eastern posi- 



tions, where the galaxy was located near a minimum of the 
sidelobes, whereas the variation was larger for the SW po- 
sition and maximum for the NW position where the max- 
imum of the sidelobe passed over the galactic disk. As a 
conclusion, the differences between the western Effelsberg 
and VLA spectra can in principle be explained by a mix- 
ture of effects (i), (ii) and (iii). However, we were not able 
to sort out which effect dominates. The north-eastern posi- 
tion only shows disk emission, whereas the south-western 
and south-eastern positions contain additional low sur- 
face density, low velocity gas (at ~ —200 kms -1 ). The 
most notable difference between the spectra occurs at the 
north-west position, where the Effelsberg and VLA spec- 
tra show three peaks: one due to galactic rotation at high 
velocities, a major peak at —50 kms -1 and a third, dis- 
tinct peak at ~ —180 kms -1 . At the latter two velocities 
gas clumps were detected by Phookun et al. (1993) in the 
north-western tail. Our spectrum reveals a bimodal ve- 
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Table 3. Parameters of the different galactic components. 



component 


N M (Mq 


) I (PC) 


halo 


16384 9.2 10 b 


1200 


bulge 


16384 3.2 10 5 


180 


disk 


16384 1.6 10 6 


240 


companion 


11000 9.2 10 6 


1200 



locity structure of the extended, low surface density gas 
tail, which is already present, though less broad, in the 
VLA synthesized spectrum. We find a total line flux of 
3.1 Jykms- 1 (2.1 10 8 M Q ) or 4.6% of the total Hi flux 
in the extended tail, more than the 3% of the total flux 
found by Phookun et al. (1993). 





Fig. 1. Hi observations of NGC 4254 (Phookun et al. 
1993). Masks were produced from the data cube which 
was smoothed to a resolution of 4'. The spectra were inte- 
grated within these masks. Top panel: Hi distribution. The 
contour levels in all images are (0.02, 0.04, 0.06, 0.08, 0.1, 
0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0) times the maximum 
of the image. Middle panel: velocity field. The contour lev- 



4. The model 

Since the model we used is described in detail in Vollmer 
et al. (2001) and Vollmer (2003), we will only briefly 
summarize its main features here. The particles represent 
gas cloud complexes which are evolving in an analytically 
given gravitational potential of the galaxy. 

The non-collisional component of the NGC 4254 
model consists of 49 125 particles, which simulate the 
galactic halo, bulge, and disk. The characteristics of the 
different galactic components shown in Tab. are the 
number of particles N, the particle mass M, and the 
smoothing length I. The particle trajectories are inte- 
grated using an adaptive timestep for each particle as de- 
scribed in Springel et al. (2001) and applying the following 
criterion for an individual timestep: 



At; 



20 kms~ 
a; 



(1) 



where a, is the acceleration of particle i. The minimum of 
all At; defines the global timestep used for the Burlisch- 
Stoer integrator that integrates the collisional component. 

The initial conditions setup was made using the pro- 
gram described in Hernquist (1993) with the following pa- 
rameters: 

— halo: mass M = 1.5 10 11 M , core radius r c —3 kpc, 
tidal radius r t =30 kpc, cutoff radius R=90 kpc, 

— bulge: mass M = 5.2 10 9 M Q , scale length /=0.3 kpc, 
cutoff radius i?=30 kpc, 

— disk: mass M = 2.6 10 10 M©, scale length 1=3 kpc, 
cutoff radius i?=45 kpc, disk thickness zo=600 pc. 

The collisional component of the model consists of 
20 000 model particles representing gas cloud complexes 
of different masses that are rotating within the gravi- 
tational potential of the galaxy. The total gas mass is 
= 5.8 10 9 Mq, which is -30% larger than the ob- 
served gas mass (Phookun et al. 1993). To each parti- 
cle a radius is attributed depending on its mass. During 
the disk evolution the particles can have inelastic colli- 
sions, the outcome of which (coalescence, mass exchange, 
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or fragmentation) is simplified following Wiegel (1994). 
This results in an effective gas viscosity in the disk. 

As the galaxy moves through the intracluster medium 
(ICM), its clouds are decelerated by ram pressure. Within 
the galaxy's inertial system its clouds are exposed to a 
wind coming from the opposite direction of the galaxy's 
motion through the ICM. The temporal ram pressure pro- 
hie has the form of a Lorentzian, which is realistic for 
galaxies on highly eccentric orbits within the Virgo clus- 
ter (Vollmer et al. 2001). The effect of ram pressure on the 
clouds is simulated by an additional force on the clouds in 
the wind direction. Only clouds which are not protected 
by other clouds against the wind are affected. 

To model NGC 4254, the 20 000 collisional particles 
with a 1/R column density profile were added to the 49 125 
non-collisional particles and the system was evolved dur- 
ing 2 Gyr in order to let it reach a relaxed state. At 
the end of the simulation the difference between the to- 
tal energy and the total angular momentum of the sys- 
tem was smaller than 0.5% of their initial values. During 
the last Gyr of this simulation the disk properties (disk 
scale height, surface density profile, and density profile) 
changed less than 5%. The final state of this simulation 
was used as the initial state for the simulations presented 
here. 

The companion galaxy with a mass of 10 11 M Q is 
simulated by 11000 non-collisional particles forming a 
Plummer sphere with a core radius of 2 kpc. 

5. The simulations 

5.1. The tidal interaction 

We have modeled a close and rapid encounter between 
NGC 4254 and another massive Virgo cluster galaxy, for 
which we will try to identify a possible candidate at the 
end of this section. The closest distance between the two 
galaxies is D = 27 kpc, which occurs at t = 300 Myr with 
a velocity difference of 710 kms^ 1 . Snapshots of this tidal 
interaction are shown in Fig. |21 The position and inclina- 
tion angle of NGC 4254 are used. The azimuthal angle is 
fixed in a way to reproduce the observed spiral structure. 
The galaxy rotates clockwise. NGC 4254 starts building 
a spiral structure ^100 Myr after the closest passage of 
the perturbing galaxy. The spiral arm at the western side, 
where the interaction happened, is more prominent than 
the one at the opposite side. After t = 500 Myr an asym- 
metric spiral pattern is formed with a prominent m=l 
mode. 

The final location of the perturbing galaxy in our 
simulation is 9 kpc east and 120 kpc north of the cen- 
ter of NGC 4254, at which point its radial velocity 
is —550 kms -1 with respect to the systemic velocity 
of NGC 4254. A possible candidate is the SB0 galaxy 
NGC 4262, which is located 51 kpc east and 137 kpc north 
of NGC 4254. Its radial velocity is ~ -1050 kms -1 with 
respect to the systemic velocity of NGC 4254. NGC 4262 
has a ring of Hi around it with twice the diameter D25 



of the optical disk and an inclination at least 50° with re- 
spect to the plane of the optical disk (Krumm et al. 1985). 
With an Hi rotation velocity of 220 kins" 1 at 2' radius its 
total mass is near the 10 11 M Q used in our model. Since 
NGC 4262 is an SO galaxy and therefore much more con- 
centrated than a spiral galaxy, we do not expect to ob- 
serve an impact of the interaction on its stellar content 
(see also NGC 4435 in the NGC 4438/NGC4435 system, 
Combes et al. 1988). Concerning the Hi distribution of 
NGC 4262, outer ring structures are frequently observed 
in gas-rich SO galaxies, see e.g. van Driel & van Woerden 
(1991), who considered that, in general, their origins are 
most likely a galaxy merger or a star formation burnt-out 
scenario. Among the SO galaxies studied by van Driel & 
van Woerden, NGC 4262 has by far the highest inclina- 
tion difference (60°) between the planes of its outer Hi ring 
and the stellar disk. For NGC 4262 we can only speculate 
that the interaction with NGC 4254 is at least partly re- 
sponsible for the peculiar ring structure. Its approximate 
symmetry is consistent with the time delay between the 
gravitational interaction and today, which is about one 
rotation period of the ring. 

5.2. Tidal interaction and ram pressure 

In a second step we introduce a temporal ram pressure 
profile of the form: 

_ 10 1 2 /^ 2 , } 

Plam _ (f - t ) 2 + t\ ' [ ' 

where po = 10 -4 cm -3 , vq = 1000 kms -1 , to = 1 Gyr, 
and t\ — 115 My. Since the inclination angle between the 
disk and the galaxy's orbital plane within the cluster is 
70°, the stripping occurs nearly face-on and is very ef- 
ficient - see Vollmer et al. 2001 for a discussion of the 
stripping efficiency with respect to the inclination angle. 
While face-on stripping (i > 60°) and thus a high stripping 
efficiency is necessary to reproduce the observed north- 
western very low surface brightness plume, the inclina- 
tion angle on the other hand must be smaller than 90°, 
i.e. face-on, to produce the north-west asymmetry of the 
plume. For these simulations the shape of the galaxy's or- 
bit within the cluster does not play a dominant role. The 
important conditions are (i) that the galaxy is entering the 
cluster and (ii) that ram pressure at the time of interest is 
Pram(580 Myr) = 0.7/9o^ 2 , consistent with its radial veloc- 
ity and the ICM density derived from X-rays (Schindlcr 
et al. 1999). 

However, when we run our simulations with a constant 
cloud column density of So = 10 21 cm~ 2 , this ram pres- 
sure does not produce a gas tail as observed for NGC 4254, 
while increasing the ram pressure to produce such a gas 
tail has too great an effect on the gas disk, and pushes 
a considerable fraction of the south-eastern part of the 
galaxy towards the north-west, which is not observed. 

Therefore, we decided to add the effects of cloud ex- 
pansion and evaporation to the numerical code: once the 
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Fig. 3. Snapshots of the evolution of the NGC 4254 system model, undergoing only a tidal interaction with another 
massive Virgo cluster galaxy. The left panels show the stellar disk, the right panels show the gas. The timesteps are 
marked on the upper right of each panel showing the gas. The center position of the perturbing galaxy is marked as 
a cross. 



clouds are pushed out of the galactic disk by ram pres- 
sure, they may expand or evaporate within the hot intra- 
cluster medium (Cowie & McKee 1977). To simulate these 
effects, the column density of those clouds that are lo- 
cated in regions with a local density below a thresh- 
old p\imi is made to decrease exponentially with time as 
£ = So exp(— t/t^j, where ti = 50 Myr. This corresponds 
approximately to the timescale of saturated evaporation 
(Vollmer et al. 2001). The threshold gas density is fixed 
at pii m = 10~ 2 cm~ 3 , and the parameters t% and p\\ m are 



chosen such that the Hi properties of NGC 4254 are re- 
produced. Snapshots of the evolution of the NGC 4254 
system undergoing a tidal interaction as well as ram pres- 
sure stripping are shown in Fig. The position and incli- 
nation angle of NGC 4254 are used. The evolution of the 
stellar content is very similar to that in the previous sim- 
ulation without ram pressure stripping (Fig. HJ, because 
ram pressure affects selectively only the gas. From the 
time of closest passage (t — 300 Myr) onwards, there is an 
increase in the difference between the gas distribution in 
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Fig. 4. Snapshots of the evolution of the NGC 4254 system model undergoing as tidal interaction with another massive 
Virgo cluster galaxy as well as ram pressure stripping. The left panels show the stellar disk, the right panels show 
the gas. The timesteps are marked on the upper right of each panel showing the gas. The center position of the 
perturbing galaxy is marked as a cross. The arrow indicates the direction of ram pressure, and the length of the arrow 
is proportional to its strength. 



the tidal interaction simulation and the mixed interaction 
simulation. Since the gas clouds that are expelled from the 
galaxy by the tidal interaction have a lower local gas den- 
sity, their column density decreases, and they are stripped 
easily by a moderate amount of ram pressure. This can be 
noted by the amount of low volume density gas that has 
been pushed towards the north-west at t — 400 Myr. In 
the forming western spiral arm the high local gas density 
prevents the gas from being stripped, and only the gas 
located outside the arm is stripped. 



5.3. Comparison between observations and simulations 



We chose the snapshot at t = 580 Myr for the com- 
parison with the observations, because (i) it reproduces 
the spiral pattern of NGC 4254 and (ii) it shows a very 
low averaged surface density tail to the north-west as it 
is observed. At that moment the ram pressure is only 
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Pram = 0. 7(0(^0 = ^0 cm 3 (km s x ) 2 . Using the Gunn & 
Gott (1972) approximation we obtain: 



£< t iT 



Pram j 



(3) 



where S is the gas surface density, v TO t is the rotation ve- 
locity of the galaxy and R is the galactic radius. With 
v mt = 170 kms -1 and £ = 10 20 cm -2 we obtain R = 
14 kpc, or about the optical size D 2 5 of the galaxy. Thus 
ram pressure affects only gas with £ < 10 20 cm -2 . Since 
the radius where the radially averaged £ = 10 20 cm~ 2 of 
the Hi gas in NGC 4254 is 18 kpc (Warmels et al. 1988, 
Cayatte et al. 1994), i.e. 1.2 times the optical radius, we 
do not expect that the galaxy would be affected by ram 
pressure. We suggest that the peculiar low surface density 
gas distribution and kinematics are due to the interplay 
between the tidal interaction and ram pressure stripping. 
The observed and the simulated distribution of the stellar 
content are shown in Fig. [3] Like the observed main spiral 
pattern, the simulated spiral structure is asymmetric with 
a wider, western spiral arm and a tighter, eastern spiral 
arm. In addition, the observed spiral structure shows an 
m=3 and an m=5 mode, which are not reproduced by the 
model. The excitation of these modes depends crucially 
on the mass distribution within the galaxy. Since the aim 
of this study is not to reproduce exactly the complicated 
spiral pattern of NGC 4254 but its gas distribution and 
kinematics, we have not attempted to investigate the influ- 
ence of different mass distributions on the spiral pattern, 
which also would require a larger number of particles. 
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5.4. The moment maps 

For the moment maps of the Hi data we smoothed the 
cube to a resolution of 35" x 35" and clipped the data 
at a 0.8 mjy/beam/channel level. A comparison between 
simulated and observed gas distributions at t = 580 Myr 
is shown in Fig. |SJ Since we do not differentiate between 
atomic and molecular hydrogen, the simulated gas distri- 
bution shows the total gas content. Although the model 
gas distributions are more concentrated in the central part 
of the galaxy than is observed in Hi, it should be noted 
that in the central few kpc of NGC 4254 most of the gas is 
molecular (Kenney & Young 1988), which makes the to- 
tal gas distribution resemble the simulations. In the sim- 
ulations without ram pressure stripping, the gas follows 
closely the stellar spiral structure, as expected. This is 
also the case in the inner parts of the galaxy when ram 
pressure stripping is included, whereas the outer parts of 
the gas distribution show important differences, since ram 
pressure has pushed the gas from the south-eastern quad- 
rant towards the center. We suggest that the observed lo- 
cal Hi maximum in the south-east is due to this effect. The 
low surface density northern extension is fairly well repro- 
duced by the simulations including ram pressure stripping. 
However, this feature extends more to the west in the ob- 
served data. The observed north-western, very low surface 
density plume is present in the simulations including ram 
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Fig. 5. The stellar component of NGC 4254. Upper panel: 
2MASS K band image (from NED). Lower panel: model 
distribution, for the model with both tidal interaction and 
ram pressure stripping; snapshot at t = 580 Myr. 



pressure stripping, but their surface density structures dif- 
fer. Whereas in the simulations there is a gradient to lower 
surface densities to the north, the observations show a sud- 
den drop between the galactic main disk and the plume 
(see also Fig. P). Within the plume the surface density 
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seems to be rather constant. The observed southern low 
surface density blob is not present in our simulations. 

A comparison between the simulated and observed ve- 
locity fields at t = 580 Myr is shown in Fig. [7J The effects 
of ram pressure stripping on the velocity field are the fol- 
lowing: (i) on the eastern side of the galaxy the gas with 
low radial velocities is pushed towards the north-west, and 
some of the gas is decelerated towards the cluster mean 
velocity (1100 kms -1 ). The net effect is that the eastern 
maximum is displaced to the south and has a more nar- 
row shape; (ii) the western plateau of the rotation curve is 
less pronounced when ram pressure stripping is added to 
the simulations; (iii) the velocity contours are straightened 
to the north, but the imprint of the spiral arm's velocity 
structure is still visible. This velocity structure is more 
pronounced in the observed velocity field. The straight- 
ening of the model velocity contours to the north can be 
observed in Fig. ^ The observed negative velocities with 
respect to the systemic velocity at RA=12 h 18 m 50 s and 
DEC=14°28'00" do not have counterparts in the simula- 
tions. 

The comparison between simulated and observed ve- 
locity dispersions at t = 580 Myr are shown in Fig. [S] The 
velocity dispersion of the simulations without ram pres- 
sure stripping show only a significant enhancement in the 
galaxy center, as expected. When ram pressure stripping 
is included, three additional local maxima are found: (i) to 
the east, (ii) to the north, and (iii) to the west of the cen- 
ter. The observed velocity dispersions also show extended 
maxima to the east, north andwest of the galaxy center 
(Fig. m . 

5.5. The three dimensional data 
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Fig. 6. The gaseous component of NGC 4254. Upper 
panel: model distribution of the Hi gas for the model with 
only a gravitational interaction. Middle panel: model dis- 
tribution of the Hi gas for the model with gravitational in- 
teraction plus ram pressure stripping. Lower panel: Hi dis- 
tribution of NGC 4254 observed by Phookun et al. (1993). 
The contour levels in all images are (0.02, 0.04, 0.06, 0.08, 
0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0) times (i) half 



Since the three dimensional Hi distribution of NGC 4254 
has a complex structure, moment maps do not show all 
the information contained in the data cubes. To complete 
the information contained in the moment maps, we show 
some selected position- velocity diagrams as well as a three 
dimensional visualisation of the model and observed Hi 
cubes. 

Fig. shows the Hi data and the model including a 
gravitational interaction and ram pressure in 3 position- 
velocity diagrams along the declination axis and 3 others 
along the right ascension axis. The Hi data is represented 
by contours and the model by the greyscale. The observed 
main features are all qualitatively reproduced. However, 
the low velocities of the low intensity, south-eastern edge 
(AS = -120", 0"; Aa = 160") are not present in the 
model. We speculate that this represents gas initially lo- 
cated at the tip of the spiral arm that is decelerated by 
ram pressure to smaller velocities. However, a captured 
dwarf scenario might also be possible, which affects the 
gas kinematics locally. No optical counterpart is detected 
at this location. In addition, the model gas distribution at 
the northern edge (AS = 200") is offset from the observed 
gas distribution as it can be seen in Fig. El 
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Fig. 9. Position-velocity diagrams of 6 selected planes. Contours: Hi observations (Phookun et al. 1993). Greyscale: 
model including ram pressure. Dark regions correspond to regions of low intensity. 



A three dimensional visualisation of the full data set 
is shown in Fig. 1101 where the upper and middle panels 
represent, respectively, the final snapshots of the simula- 
tions (i) with a gravitational interaction alone (cf. Fig. [3J| 
and (ii) with a gravitational interaction and ram pres- 
sure (cf .Fig. |3J|. The lower panel shows the Hi data 
cube of Phookun et al. (1993) with an isosurface at 
2 mjy/beam/channel. The two tidally induced spiral arms 
are clearly visible to the left and right of the upper panel. 
By comparing the upper and lower panel the effects of 
ram pressure on the model gas distribution can be seen di- 
rectly: the extent of the gas distribution in right ascension 
is decreased (cf. Fig.EJ, the velocities in the high velocity 
arm are increased (the gas is accelerated; the wind and 
galactic rotation have the same direction) and those in 
the low velocity arm are decreased (the gas is decelerated; 
the wind and galactic rotation have opposite directions). 
These features can be qualitatively recognised in the Hi 
data cube fFig. ITU1 lower panel). In addition, the quasi 
horizontal edge at v ~ 2460 kms -1 to the right of Fig. 1101 
(upper panel) is stretched to higher velocities in Fig. 
(middle panel) due to the gas acceleration by ram pres- 
sure. However, we note that the observed gas velocities 
in Fig. ^3 (lower panel) at that location in the data cube 
are shifted to much higher velocities (Av ~ 100 kms -1 ). 
Moreover, the low velocities in the Hi blob located in the 
upper left quadrant of Fig. HTH flower panel), which corre- 



sponds to the southern blob of Fig.EI are not reproduced. 
If we lower the intensity level of Fig. (middle panel) , 
single gas clouds become visible in this region of the cube, 
but they do not form an elongated structure as it is ob- 
served. 

Thus we conclude that although the observed main 
features are reproduced qualitatively, there are quanti- 
tative differences between observations and simulations. 
However, we note that effects of ram pressure definitely 
improve the fit to the Hi observations (gas displacement, 
acceleration and deceleration) . Only the amplitude of the 
mechanism is not strong enough to reproduce the Hi in 
the lower right and upper left quadrant of Fig. 1101 (lower 
panel). We speculate that this is mainly due to the miss- 
ing physics of Hi cloud expansion and evaporation in the 
model and, in a less important way, to the amount of ram 
pressure and the inclination angle between the disk and 
the galaxy's orbital plane within the cluster. 

5.6. The spectra 

The observed and simulated Hi spectra of the main body 
of the Hi distribution (see Sect. 13 are shown in Fig. ^] 
While all spectra show the characteristic double-horned 
structure, only the simulation that includes ram pressure 
stripping reproduces the observed tail at low velocities. 
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Fig. 7. Top panel: model distribution of the Hi veloc- 
ity field for a model of a gravitational interaction alone. 
Middle panel: distribution of the Hi velocity field of the 
model of a gravitational interaction together with ram 
pressure. Bottom panel: Hi velocity field of NGC 4254 
(Phookun et al. 1993). The contour levels of all images 



Fig. 8. The gaseous velocity dispersion in NGC4254. 
Upper panel: model gas velocity dispersion map for the 
model with only a gravitational interaction. Middle panel: 
model gas velocity dispersion map for the model with 
both gravitational interaction and ram pressure stripping. 
Lower panel: Hi velocity dispersion map (Phookun et al. 
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Fig. 10. 3D representation of data cubes. The axis of right 
ascension is indicated. The axis of declination is perpen- 
dicular to the image plane. Upper panel: model cube of 
the simulation including a gravitational interaction alone. 
Middle panel: model cube of the simulations including a 
gravitational interaction and ram pressure. Lower panel: 
observed Hi data cube (Phookun et al. 1993). The inten- 



Fig. 11. The central Hi spectrum of NGC 4254. Upper 
panel: Simulated central Hi spectrum at t = 580 Myr for 
the model with only a gravitational interaction. Middle 
panel: Simulated central Hi spectrum at t = 580 Myr of 
the model with both gravitational interaction and ram 
pressure stripping. Lower panel: Observations. Solid lines: 
Effelsberg 100-m single dish spectra of the four off-center 
positions. Their locations with respect to the galaxy center 
are marked on top of each panel. Dashed line: synthesized 
VLA spectrum. Dotted line: 3<r noise level of the 100- 
m spectrum. Radial velocities are given relative to the 
systemic velocity of NGC 4254. 



The observed and simulated Hi spectra of the 4 off- 
center positions (see Sect.[3J) are shown in Fig. ^] for the 
model with only gravitational interaction and in Fig. 1131 
for the model with both gravitational interaction and ram 
pressure stripping. 
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Fig. 12. Off-center positions Hi spectra of NGC 4254. 
Solid line: simulated off-center spectra at t = 580 Myr for 
the model with only a gravitational interaction. Dotted 
line: observed Effelsberg off-center spectra. The locations 
of the spectra with respect to the galaxy center (see 
Sect. are marked on top of each panel. Radial veloc- 
ities are relative to the systemic velocity of NGC 4254. 
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Fig. 13. Off-center positions Hi spectra of NGC 4254. 
Solid line: simulated off-center spectra at t = 580 Myr 
for the model with both gravitational interaction and ram 
pressure stripping. Dotted line: observed Effelsberg off- 
center spectra. The locations of the spectra with respect 
to the galaxy center (see Sect. 01 are marked on top of 
each panel. Radial velocities are relative to the systemic 
velocity of NGC 4254. 



In the model without ram pressure stripping the spec- 
trum at the NE position can already be reproduced well. 
The shift of the observed Hi spectrum to higher velocities 
at v = 100 — 150 kms -1 is due to the fact that the model 
does not reproduce the high velocities of the northern Hi 
arm (see Sect. 15 .5 [I . Although the shape of the SE spec- 



trum is reproduced, the observed flux densities are 1/3 
higher than the simulated ones. The simulated SW spec- 
trum lacks the observed emission at low velocities, and 
the observed and simulated flux densities differ by a fac- 
tor of two. The simulated NW spectrum completely lacks 
the observed peak at low velocities. 

In the model including ram pressure stripping the sim- 
ulated NE spectrum is almost unchanged and reproduces 
well the observations. The simulated SE and SW spectra 
have lower flux densities than those for the simulations 
without ram pressure (Fig. I12f> and give better fits to the 
observed spectra. The effects of ram pressure stripping on 
the NW spectrum are twofold: (i) the high-velocity max- 
imum is decreased and (ii) a low velocity tail appears - 
both effects are observed. Thus all 5 observed spectra are 
well, or, in some cases better reproduced by the simula- 
tions including ram pressure stripping. 

6. Discussion 

Previously, two explanations were proposed for the m=l 
mode of NGC 4254: (i) its intrinsic mass distribution that 
favors the formation of this spiral pattern (lye et al. 1982), 
and (ii) the excitation of the m=l mode by infalling Hi 
gas together with a swing amplifier (Phookun et al. 1993). 
Based on our simulations we suggest another mechanism: 
a close and rapid tidal interaction that triggers an m=l 
mode. Our models indicate that NGC 4254, when entering 
the Virgo cluster, had a close and rapid encounter with 
a 10 11 M galaxy ~ 280 Myr ago, and has since then 
been undergoing increasing ram pressure stripping. The 
tidal interaction caused the m=l spiral structure, and, 
together, the tidal interaction and ram pressure caused 
the peculiar Hi distribution and kinematics. The intrinsic 
mass distribution of NGC 4254 then favors the formation 
of additional m=3 and m=5 spiral structures. 

Phookun et al. (1993) discarded the ram pressure sce- 
nario on the following grounds: since NGC 4254 is mov- 
ing towards the south-east, which makes the ram pres- 
sure push the gas located in the south-eastern part of the 
galaxy towards the north-west, there can be no gas ahead 
of the galaxy as it moves into the intracluster medium. 
However, the observed low surface density gas blobs to 
the south of the main Hi body represent such gas. 

Although it is true that in an edge-on stripping sce- 
nario the existence of the southern, low surface density 
gas blobs would rule out the ram pressure scenario, in our 
model stripping occurs almost face-on. In the case of a 
face-on stripping the gas is removed in a cone-like struc- 
ture and the projected radial distance of the stripped gas 
increases with its height from the disk (see Vollmer & 
Huchtmeier 2003). This makes it in principle possible to 
observe low surface density gas seen in projection to the 
south of the main gas disk, which represents stripped gas 
that is actually located in front of NGC 4254 and which 
should show negative velocities with respect to the sys- 
temic velocity of NGC 4254 as is observed. The southern, 
low surface density gas blobs are thus not contradictory to 
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a face-on ram pressure stripping scenario, and we in fact 
find gas clouds in the south of the gas disk in the model 
snapshot that we compare to the observations in Fig.Q]- 
these are not visible in the gas distribution (Fig. [fJJ and 
3D visualisation (Fig. ITU|l. because their column density 
is very low. 

Although our model including ram pressure stripping 
succeeds in reproducing the spiral structure and the over- 
all gas distribution and kinematics, the main difference 
between our model and the observations is the shape of 
the extended, low surface density gas tail observed to the 
north-west. Whereas the simulated gas tail is a smooth 
continuation of the m=l spiral arm, the observed tail con- 
sists of single blobs, at the level of sensitivity of the VLA 
data, and there is a gap between the gas disk and the 
tail, which might have two reasons: (i) Once the stripped 
gas leaves the disk it is ionized and expanding through 
ICM-ISM mixing, and is thus no longer visible in the Hi 
line. At a later stage it might cool and recombine again, 
(ii) The stripping is not continuous, but occurs in discrete 
events, because the ICM is not homogeneous. In this case 
the observed tail might represent gas falling back. We did 
not attempt to simulate this effect, because we do not 
have enough information on the temporal profile of such 
discrete events. 

The evolutionary history of NGC 4254 thus closely re- 
sembles that of NGC 4654 (Vollmer 2003). Both galax- 
ies had a close and rapid encounter with another mas- 
sive Virgo cluster galaxy and both galaxies are presently 
still undergoing ram pressure stripping. In both cases, the 
tidal interactions make the stripping more efficient. For 
NGC 4254 we had to include the effects of gas expansion 
and/or evaporation into the simulations, which was not 
necessary for NGC 4654 - if we had included these effects 
in the NGC 4654 simulations, we could have used a ram 
pressure that increases with time and whose maximum 
value would have been a few times smaller. 

Recently Minchin et al. (2005) detected 1.6 x 10 8 M 
of neutral hydrogen in the object VIRGOHI21 which has a 
radial velocity of ~ 2000 kms -1 and lies ~ 120 kpc north- 
west from NGC 4254. The total Hi flux of VIRGOHI21 is 
too low and its position too far away from the beam center 
of our NW position to be detectable in the Effelberg spec- 
trum. We exclude that VIRGOHI21 represents tidal debris 
from NGC 4254 for two reasons: (i) the velocity gradient 
within VIRGOHI21 is opposite to the gradient between 
NGC 4254 and VIRGOHI21 and (ii) if VIRGOHI21 were 
tidal debris, the second interacting galaxy should be lo- 
cated roughly on an extrapolated line between NGC 4254 
and VIRGOHI21, but there is no galaxy close enough in 
velocity in this region of the sky. 

7. Conclusions 

NGC 4254 is a peculiar spiral galaxy located at the pe- 
riphery of the Virgo cluster. It has a prominent m=l spi- 
ral pattern, making it look basically one-armed, and its 
Hi distribution and kinematics are quite distorted. Deep 



VLA Hi imaging (Phookun et al. 1993) showed kinemati- 
cally perturbed, low surface density gas around the galaxy, 
in particular, an extended, very low surface density tail to 
the north-west of the galaxy. We observed NGC 4254 in 
Hi with the Effclsberg 100-m telescope at 5 different po- 
sitions: towards the center and at 4 positions located one 
9.3' HPBW to the NE, NW, SW, and SE from the center. 
Our NW spectrum shows ~ 30% more low surface density 
gas than the VLA data, and a bimodal velocity distribu- 
tion in the part that corresponds to the north-western gas 
tail. 

We investigated if a scenario involving a tidal interac- 
tion and subsequent ram pressure stripping could explain 
these properties. 

In a first step we modeled a close and rapid encounter 
between NGC 4254 and another massive galaxy. In this 
way we were able to reproduce the observed asymmetric 
spiral structure ~ 280 Myr after the passage of the per- 
turbing galaxy, for which we identified NGC 4262 as a 
potential candidate. 

In a second step we included the effects of ram pressure 
stripping that slowly increases with time since we assume 
that NGC 4254 is entering the Virgo cluster, in an almost 
face-on stripping geometry, with an angle of 70° between 
the galaxy disk and the orbital plane. At its location at 
the periphery of the Virgo cluster, the small ram pressure 
(Pram ~ 100 cm -3 (kms -1 ) 2 ) would not cause significant 
Hi stripping of NGC 4254 without the tidal interaction 
event that has pulled out gas from the disk. The observed 
Hi properties can only be reproduced if we include the ef- 
fects of gas expansion and/or evaporation into our simula- 
tions. Combining all effects we succeeded in reproducing 
the main characteristics of the Hi spectra, gas distribu- 
tion, and velocity field, as well as the m=l spiral pattern. 

Thus, NGC 4254 had an evolutionary history in the 
Virgo cluster similar to that of NGC 4654 (Vollmer 2003); 
both are entering the Virgo cluster for the first time and 
experienced a tidal interaction and still undergo ram pres- 
sure stripping. 
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